Self-localization of polariton condensates in periodic potentials 
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We predict the existence of novel spatially localized states of exciton-polariton Bose-Einstein 
condensates in semiconductor microcavities with fabricated periodic in-plane potentials. Our theory 
shows that, under the condition of continuous off-resonant pumping and losses associated with 
polariton decay, localization is observed for a wide range of optical pump parameters due to effective 
potentials self-induced by the polariton flows in the spatially periodic system. We reveal that the 
localization of exciton-polaritons in the lattice may occur both in the gaps and bands of the single- 
particle linear spectrum, and is dominated by the effects of gain and dissipation rather than the 
applied potential, in sharp contrast to the conservative condensates of ultracold alkali atoms. 



Introduction. - Ultracold atomic gases loaded into pe- 
riodic potentials created by optical or magnetic lattices 
continue to provide a fascinating tool for exploring funda- 
mental quantum many-body physics, as well as offering 
an opportunity to build quantum simulators and scalable 
quantum information processing systems Q . 

In recent years, polariton condensates created in semi- 
conductor microcavities [2^ have emerged as an attractive 
alternative to the atomic systems, owing to the relatively 
high condensation temperatures, direct momentum and 
real-space imaging via the cavity photoluminescence, and 
possibility to manipulate a condensate using both optical 
pump and structure of the microcavity [3] . Considerable 
effort has been directed towards achieving the polariton 
condensation in one- and two-dimensional 'lattice' poten- 
tials and observing macroscopic population of the single- 
particle band-gap spectrum. So far, the majority of the 
experiments with exciton-polaritons in a periodic poten- 
tial have been performed in a pulsed excitation regime, 
which has allowed to populate and image rapidly decay- 
ing high momentum states in a periodic potential. The 
staggered V-state' phase structure of the higher-order 
Bloch states has been demonstrated in a ID lattice cre- 
ated by metal deposition onto the microcavity Sim- 
ilarly, condensation into second band has been achieved 
in a 2D square lattice [5], and the Brillouin zone struc- 
ture of the Kagome lattice has been mapped out [6] . Re- 
cently, condensation into an energy state corresponding 
to a single-particle spectral gap has been demonstrated 
in a ID lattice created by lateral modulation of quantum 
wires in a microcavity . Apart from the fabricated lat- 
tices, condensation in the tuneable acoustically created 
lattices has been investigated [H, . 

One of the most fundamental effects due to the band- 
gap structure is localization of the condensate in real 
space. In Kagome lattices [6] this, in principle, can be 
achieved by accessing the flat band of the spectrum, 
which signifies dispersionless dynamics. Alternatively, 
spatial localization in spectral gaps of the band-gap 
structure can be achieved due to the balance of repul- 
sive nonlinearity and modified effective dispersion due to 



the negative effective mass. This mechanism has been 
explored in the context of gap solitons both in photonic 
[10] and matter-wave [llj band-gap structures. Scarce 
theoretical descriptions of exciton-polaritons in periodic 
potentials [H, [l3| so far disregarded the fact that, un- 
like the atomic condensate in thermal equilibrium, the 
polariton condensate is a driven, highly dissipative sys- 
tem due to finite lifetime of quasiparticles and the pres- 
ence of the optical pump. This distinct dissipative nature 
of polariton condensates becomes very important in the 
continuous- wave (cw) non-resonant excitation regime (H , 
which creates long-living condensates with spontaneously 
established spatial coherence and long coherence times. 
The presence of gain-assisted localization and related dy- 
namics may play a critical role for such dissipative con- 
densates [lJ,ll5|, in contrast to atomic Bose-Einstein con- 
densates (BECs) in thermal equilibrium. 

In this work, we identify and analyse the mechanism 
for localization of a dissipative polariton condensate in 
a non-resonant cw excitation regime, with an external 
periodic potential. We show that the driven-dissipative 
nature of the polariton condensate gives rise to addi- 
tional effective potentials local to the condensate, namely 
a defect-type potential introduced by the pump-induced 
polariton reservoir and a potential self-induced by the 
strongly modulated flow of polaritons away from the ex- 
citation spot. Trapping of the polariton condensate in 
these effective potentials, which supersede the external 
lattice potential^ is responsible for its spatial localization. 

Model. - We assume that the near-equilibrium polari- 
ton condensate is formed via cw incoherent excitation 
with a spatially inhomogeneous pump. Above conden- 
sation threshold, it can be described by the mean field 
open-dissipative model for a polariton condensate cou- 
pled to an incoherent polariton reservoir density [16]: 
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Vnl + ARuR - 7c) + Vl 



-(7fi + J?|M/|Vii(r,i) + P(r). 
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Here ^ is the condensate wavefunction, ^ ^i^" 
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FIG. 1: Steady state energies of the polariton BEG obtained 
with cw on-site (closed circles) or 7r/4 off-site (open circles) 
narrow Gaussian pump with the width cr^ = 0.6. Shaded 
(open) areas show bands (gaps) of the single-particle Bloch- 
wave spectrum in the ID lattice potential Vl{x) (Vb = 5.0, 
ki — 1.5). Dashed lines are linear defect states in the com- 
bined effective potential Vl{x) + gRP{x)/jR (see text). 



ear periodic potential created for the condensate, and 
Vnl = ^d^P + gR^R is an effective nonlinear poten- 
tial depending both on the condensate density and the 
reservoir density, ur^ which is determined by the inho- 
mogeneous optical pump P(r). The parameters defin- 
ing the condensate dynamics are the relaxation rates 
of the condensed polaritons jc and reservoir polaritons 
7i?, the stimulated scattering rate i?, and the polariton- 
polariton interaction strengths gR^ We have writ- 
ten the model ([T]) in a dimensionless form by using the 
characteristic scaling units of time l/7c, energy fi7c, and 
length [h/ {mLPlc)Y^'^ ^ where m^p is the lower-polariton 
effective mass [17]. Under this scahng 7^ = 1, however 
we formally retain this parameter hereafter. The scaling 
units are calculated using the parameters close to those 
of the experimental setup of ^Jl], with ttilp = 10~^me, 

= 6 X IQ-^ meViim^, gR 2gc, 7c 0.02 
^R = 1.57c, and R = 0.01 /im^ ps~^. 

localization in ID. - We begin our consideration with 
a simplest, one-dimensional (ID) model of a quasi-lD 
polariton condensate (e.g., confined in a microwire [l9| ) 
subject to a periodic potential approximated by a har- 
monic function Vl{x) = Vb sin^(/c^x), with the lattice con- 
stant a = ir/ki. The periodicity of the potential imposes 
a band-gap structure onto the LP branch of the polari- 
ton dispersion curve [3], and we can consider population 
of this band-gap structure by polaritons depending on 
the pump intensity and width. To this end, we look for 
the steady state solutions ^(x) = iIj{x) ex.-p{iEt) of the 
model ^ governed by the stationary complex Ginzburg- 
Landau equation with inhomogeneous nonlinear gain and 
linear loss. We will consider the inhomogeneous Gaus- 
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FIG. 2: Spatial profiles of the condensate density ln(y9(x)) 
and pump intensity ln(P(x)) in a periodic ID lattice potential 
Vl(x) with Vb = 5.0 and ki = 1.5, corresponding to the peak 
pump intensity Po/Pth = 1-27 [point (a) in Fig. 1], Po/Pth = 
2.53 [point (b)], Po/Pth = 5.06 [point (c)] , and Po/Pth = 4.30 
[point (d)]. In addition, panel (a) shows the lattice potential 
In(VL) (note the difference in the position scale). 



sian pump profile P{x) = Pq exp(— a^x^) and relate the 
peak intensity of the pump to the threshold homogeneous 
pump intensity Pth = JrJc/R [16]. 

As has been established both experimentally (see, e.g., 
and theoretically [H, HqI, cw non-resonant 
pump leads to spatial localization of polariton BEC even 
without a trapping potential. This effect, resulting in 
formation of dissipative solitons [Ti, [2l| is due to the 
balance of gain and loss in the system and is irrespective 
of the sign of the nonlinear interactions. This mecha- 
nism is dramatically different from conventional balance 
of nonlinearity and dispersion (diffraction) responsible 
for formation of bright solitons in the coherent fields. 
The repulsive polariton-polariton interaction then leads 
to considerable spatial extension of the polariton con- 

i3,[i 



,[13, [Til, and 



densate beyond the excitation spot 
results in density deformations which are esMcially pro- 
nounced in the presence of a trap [lH, [13, 112] • In what 
follows, we will demonstrate that spatial localization of 
the condensate due to the balance between gain and loss 
persists in the presence of a periodic potential. 

Narrow excitation. - A spatially narrow excitation, 
P(x), can selectively target a single lattice site and there- 
fore we can distinguish between an on-site and off-site 
excitation. In the later case, the peak of the excitation 
spot can be arbitrarily offset from the minimum of the 
lattice potential. As can be seen from examining Eq. ([1]), 
in the case of the pump intensity near the threshold for 
the condensate formation, Po/Pth ^ 1, nonlinear effects 
can be neglected. In this case, assuming a straightfor- 
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FIG. 3: Spatial profiles of the lattice potential Vl{x) (Vb = 
5.0, ki = 1.5) and the effective linear potential due to the 
phase gradient /2 corresponding to the peak on-site pump 
intensities (a) Po/Pth = 1.27, (b) Po/Pth = 5.06, and (c) a 
non-centered pump (7r/4 off-site) with Po/Pth = 4.30. Panel 
(d) shows the /c-space spectrum of the steady state [(a) in Fig. 
[2] corresponding to the effective potential in (a). 



ward balance of gain and loss terms, Eq. ([T]) would de- 
scribe a condensate in a linear periodic potential Vl{x) 
with the addition of a linear defect induced by the pump, 
gRn% ^ {gR/-fR)P{x), where n% = P/{-fR^R\ip\^) is the 
steady state reservoir density. The energy of the lowest 
order linear defect modes is plotted as a function of pump 
intensity in Fig. [H together with the band-gap spectrum 
imposed by the linear periodic potential. 

Taking this spectrum as a frame of reference, we can 
examine the energy and spatial structure of the localized 
steady state depending on the pump intensity. The re- 
spective dependence is plotted in Fig. [T]and shows that, 
near the threshold on-site excitation, the energy of the 
localized state is indeed very close to that of a symmet- 
ric linear defect mode. As the ratio Po/Pth grows, so 
does the energy of the localized state until it crosses into 
the band of the linear band-gap spectrum. This is in 
contrast to dissipative photonic lattice solitons ^3| and 
defect states [24] , which have been previously found to ex- 
ist exclusively in the gaps of the band-gap spectrum. An 
additional feature of the polariton steady states is tran- 
sition, with growing pump intensity, from narrow states 
localized in a vicinity of a single potential well [Fig. Efa)] 
to broader states with exponentially decaying and spa- 
tially modulated tails [Fig. [2fb,c)]. 

To understand the spatial structure of the localized 
states, we perform the standard Madelung transforma- 
tion iIj{x) = p{x) exp[i^(x)], and write down the steady 
state equations for the condensate density $ = ^/p and 



phase gradient (flow velocity) v = (j)^'- 

-2(vl^ 9c^^ + QRul + ^ + 2^^ = 0, (2) 

{^\), - {Rn% - 7c)^' = 0, 

Given the rapid decay of the inhomogeneous pump at 
large x, we can immediately see that the asymptotic be- 
haviour of the condensate near-linear density away from 
the pump spot is determined by the combination of the 
linear periodic potential Vl{x) and a self-induced effec- 
tive potential due to the phase gradient, v'^{x)/2. These 
two potentials for the narrow and wide states (a) and 
(b,c) in Figs. I1|2I are shown in Fig. [3fa-c). The strong 
contribution of the Bloch states from the edges of the 
Brillouin zone into the composition of the localized state 
[Fig. Efd)] leads to steep spatial gradients of the phase, 
0(x). As a result, the flow of the condensate out of the 
excitation spot induces a periodic, tight-binding Kronig- 
Penney-like potential [Fig. [3] (a)], and the band-gap spec- 
trum due to the linear lattice potential becomes of little 
relevance to the energy eigenstates. The sharp peaks 
of the effective potential are offset from the maxima of 
the lattice potential, which leads to the fact that the 
condensate density peaks can coincide with the maxima 
of the lattice potential [see Fig. Efa)]. The tunnelling 
of the condensate out of the excitation spot is strongly 
inhibited by the tight-binding effective potential, until 
the modulations of the flow are smoothed out by the 
stronger pump. As demonstrated in Figs. [2fb,c) the 
strongly pumped condensate density then extends over 
many wells of the periodic potential. We stress that, in 
contrast to the lattice potential which exists regardless of 
the condensate formation, the self-induced potential can 
become physically meaningless far on the tails of the con- 
densate [shaded areas in Fig. [3] (a)], where the density 
is negligible and any dynamical fluctuations can quickly 
destroy regular phase structure. 

Although the narrow gap states [Fig. [2fa)] greatly re- 
semble well localized gap solitons of a conservative BEG 
[111] , they have no connection to the Bloch states at the 
band edge of the first Brillouin zone in the linear lattice 
potential, rather branching off the linear defect states. 
Their spatial extent is determined purely by the phase 
gradients established due to the balance of gain and loss 
in the system. The energy of the localized state can lie 
in either a spectral gap or band of the lattice potential 
since it is supported by the effective potentials created 
by the polariton flow. 

In the off-site excitation regime, the localized states 
can also be formed both within gaps and bands of the 
linear band-gap spectrum of the potential Vl{x) [Fig. 
[1], and the threshold for the formation of these states 
is higher than that of the on-site state. Remarkably, 
the asymmetry of the defect created by the offset pump 
translates to the asymmetry of both the effective poten- 
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tial generated by the polariton flow and the condensate 
density [Fig. [2fd) and Fig. [3fc)]. Such asymmetric 
locahzed states are supported by the interface between 
two different effective potentials, and therefore resemble 
'mixed-gap' interface states of an atomic BEC in anhar- 
monic lattices (iij . 

Wide excitation. - By flxing the pump intensity and 
extending the width of the pump spot signiflcantly be- 
yond the single lattice well, we flnd a variety of spa- 
tially extended nearly energy-degenerate localized states. 
One example of such a gap state and its supporting ef- 
fective potential is presented in Figs. |4l^a,b). These 
states strongly resemble 'flat gap states' formed by the 
nonlinearity-induced truncation of nonlinear Bloch states 
in a conservative condensate [26|. However, the localiza- 
tion in this case is due to the effective potential created 
by the polariton flow out of the flnite gain region. 

A very wide, locally homogeneous pump near conden- 
sation threshold generates extended states that locally 
resemble a ground state Bloch wave [Fig. |4] (c)] of the 
linear lattice potential Vl{x). In this regime, the asymp- 
totic approximation of the phase proflle at x ^ 0, derived 
in [17} for the case of untrapped polariton condensate be- 
comes valid for the most of the excitation region [Fig. |4] 
(c), dashed line], and the effective potential due to the 
polariton flow becomes harmonic: 
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The density proflle is a typical ground state supported 
by the combination of the periodic and harmonic effec- 
tive potential [Fig. |4](c)], and its /c-space spectrum [Fig. 
|4](d)] displays localization around k = characteristic of 
the ground state in the lattice [4]. Thus, in this regime 
the dissipative condensate in a periodic potential mimics 
the behaviour of a conservative condensate in the com- 
bined periodic and harmonic potentials. 

Localization in 2D. - Similarly to the case of ID pe- 
riodic potential, we flnd that excitation spots of varying 
width can create narrow and wide localized steady states 
of the polariton condensate both in the gaps and bands 
of the single-particle band-gap spectrum corresponding 
to a 2D 'square' lattice potential Vl = Vb[sin^(/c^x) + 
sin^(/c^?/)]. Examples of a narrow state density and its 
phase are shown in Fig. [5] for the case of a weak lattice. 
Although the detailed description of such states is beyond 
the scope of this work, we note that the phase structure 
displays sharp domain walls associated with strong phase 
gradients, similar to those in ID case, which leads to the 
tight spatial localization of the condensate. 

Conclusions. - We have demonstrated that, under the 
conditions of spatially inhomogeneous non-resonant cw 
excitation, a polariton condensate can be spatially local- 
ized in a periodic in-plane potential owing to the polari- 
ton flows established due to the balance between gain 
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FIG. 4: (a,c) Spatial profiles of the condensate density 
phase (j){x), and lattice potential Vl{x) (Vb = 5.0, ki = 1.5) 
corresponding to the peak on-site pump intensity Po/Pth = 
1.27 and widths = 6x lO"'^ and = 6x 10~^, respectively. 
Dashed curve in (c) is the analytical approximation of 
Eq. Q. (b) Lattice potential Vl and the self-induced poten- 
tial, v'^ corresponding to (a), (d) Spectrum of the state (c) 
within the first Brillouin zone of the lattice potential. 




FIG. 5: Spatial profiles and cross-sections of the (a,b) con- 
densate density and (c,d) phase, corresponding to a narrow 
in-gap state in a 2D lattice (Vb = 2.0, ki = 7v/4). 



and loss. Periodic potential strongly modulates the ve- 
locity of the polariton flow out of the pumping region, 
and the polaritons are self-trapped in the resulting ef- 
fective potential. This mechanism of localization ensures 
that, depending on the pump intensity, condensate can 
be formed and localized both in gaps and bands of the 
linear single-particle band-gap spectrum of the periodic 
potential. Localization of the dissipative condensate is 
therefore irrelevant to the requirement of the negative 
effective mass which is essential for the gap soliton for- 
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mation in conservative (atomic) BECs with repulsive in- 
terparticle interactions [llj . 

The predicted features of the real-space localization 
of a dissipative polariton condensate subject to a pe- 
riodic potential have strong consequences for any pro- 
posed polariton-based devices incorporating modulated 
potentials, and should be investigated further, both the- 
oretically and experimentaly. Polariton condensates in 
'lattice' traps, created in the non-resonant cw excitation 
regime explored here, are ideally suited to such studies. 
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